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ABSTFiACT 


Solid state microwave components are all nonlinear to some degree. 
The circuit analysis of these components requires a nonlinear device 
model and analytic means to extract the effect of device-circuit 
interactions. The two main techniques used for the nonlinear circuit 
analysis are frequency-domain and time-domain techniques. The 
frequency-domain simulation can be applied to circuits having weak 
nonlinearity and also these techniques takes excessive amount of time 
if the circuit has a large number of harmonics. 

This thesis work presents a time-domain technique for solving a 
class of strongly nonlinear microwave circuits or input signals having 
a large number of harmonics. The entire circuit is discretized in time 
domain, by Backward-Euler numerical technique. Equivalent circuit 
models have been used for capacitor and inductor in each time step 
resulting in a set of nonlinear equations. In the standard methods 
available for circuit analysis such as SPICE, these equations are solved 
iteratively. At each time step, these equations are linearized using 
Newton-Raphson method and solved by matrix inversion. This leads to a 
large number of matrix inversions and linearizations of the equations, 
one for each iteration. 

In the present analysis, a procedure has been developed which 
reduces the computational time significantly by solving the nonlinear 
equations directly, using a new iteration procedure. An assumed set of 
initial node voltages are iteratively updated according to an error 
criteria, based on the difference between the assumed and the 
calculated node voltages. 

This procedure involves only one matrix inversion per time step. 
However, the convergence of the iterartive process requires that the 
sum of the absolute values of the elements of any row or column 
vectors of the system Jacobian matrix be less than or equal to one. 
This condition can be satisfied by choosing a proper time step. In the 
present study proper time step for convergence has been arrived at 





experimentally . 

The method has been applied to a harmonic converter circuit. 
Equivalent circuit model for Schottky diode, step recovery diode and 
lossy dispersive transmission line also have been discussed in this 
thesis work. 
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IJNTRODUCTION 


1.1 Introduction 

The trend toward monolithic integration of microwave circuits is 
intensifying interest in computer-aided design. The interest is 
predominantly in the analysis and design of microwave circuits 
containing nonlinear devices. The optimum design of microwave circuits 
containing nonlinear solid state devices requires an accurate technique 
for predicting their performance. The techniques used for the 
nonlinear circuit analysis can be classified mainly into two categories, 
viz., frquency-domain and time-domain techniques. 

i2 Freouency-domain methods 

In frequency-domain techniques, as the name itself suggests, the 
solution to the circuit problem is effected in the Fourier transform 
domain. In the transformed domain the signal is often represented by its 
Fourier series with finite number of terms. The time-domain differential 
equations governing the circuit are converted to nonlinear algebraic 
equations in the transformed domain using Fourier techniques. These 
equations are solved for different harmonics using iterative 
procedures. 

The frequency-domain techniques are well suited for linear 
circuits, especially the distributed components. It is usually possible to 
find closed form expressions for distributed components in frequency- 
domain, therefore simulating these components is inexpensive. The 
active devices are generally nonlinear, and modelling of these 
components in frequency-domain is difficult. Also these techniques are 
inefficient if the circuit generates a large number of harmonics or the 
input has a large number of harmonics. 

E.M. Bailey [11 has introduced a frequency-domain technique for 
nonlinear system, popularly known as Harmonic balance method. In this 



method, each current or voltage variable is represented by a Fourier 
series that satisfies the requirement of periodicity. An optimization 
algorithm is then used to adjust the coefficients of the Fourier series 
such that the system equations are satisfied with least error. Its main 
disadvantage is the large number of variables that must be optimized. 
Nakhla and Vlach C23 have proposed a method which can reduce the 
number of variables to be optimized, but the method is efficient only 
for circuits having few harmonics or weak nonlinearity. 

Chua and Ushida C31 have developed algorithm for computing almost 
periodic steady state response of nonlinear systems with multiple 
frequencies. Hicks and Khan [43 have introduced a general method for 
the analysis of microwave circuits which contains a sinusoidal source 
and one or more nonlinear devices. They have achieved good 
convergence with this method. This technique cannot be applied to 
circuits having many input harmonics. 

Gilmore [53 has modified the harmonic balance technique by 
analyzing the linear elements in the frequency-domain and the nonlinear 
elements in the time-domain, the conversion between the frequency- 
domain solution of the linear network and the time-domain solution of 
the nonlinearity is usually accomplished using Fast Fourier Transform 
techniques. This limits this method to systems having only related 
signal components. 

Lamnabhi [63 has applied Volterra series technique to nonlinear 
systems. The Vol terra series relates the output to the input by a 
series of convolutions in the time-domain. Unfortunately, Vol terra 
series techniques are restricted to weakly nonlinear systems because 
of the algebraic complexity of determining Volterra nonlinear transfer 
functions of high order. 

Steer and Khan [73 used a generalized power series expansion of 
the input-output characteristic of a nonlinear system. This method can 
be used for multifrequency input systems, but requires modelling of 
devices in the generalized power series expansion and involves the 



evaluation of Jacobian matrix for large number of harmonics. 

In general, the frequency-domain techniques are inefficient if 
the nonlinearity of the circuit is high and/ or input has a very large 
number of harmonics. 

13 Time-domaw methods 

The time-domain methods find the voltages and currents in 
discrete time steps. The system differential equations are converted to 
algebraic nonlinear equations by numerical technique. The general 
purpose circuit analysis program like SPICE [83 involves these 
procedures. The characterization of a nonlinear device by a time-domain 
model is usually appropriate because circuit models typically relate 
output parameters to input ones through device physics. The distributed 
elements are difficult to handle in time-domain, since they are 
described by partial differential equations. 

Silverberg and Wing [93 have proposed a time-domain method for 
analyzing distributed parameter network containing lumped nonlinear 
elements. In this method, the linear part of the circuit is described in 
the frequency-domain and then the description is converted into the 
time-domain by numerical transform techniques. The solution of the 
whole network is obtained step by step in time at the interface of the 
two parts by solving simultaneously the convolution equations for the 
linear part and the differential equations for the nonlinear part. 
Because of the inverse-transform calculations, the method is limited to 
small sized circuits. 

1.4 Present work 

The present work on analysis of microwave nonlinear circuits 
uses time-domain method for solving a class of microwave circuits 
having a large number of input harmonics (eg. triangular pulse or step). 
In SPICE, the nonlinear equations are solved iteratively. The nonlinear 
equations are linearized using Newton-Raphson numerical technique and 



then the equations are solved using matrix inversion. This leads to a 
large number of matrix inversions per time step. In the present method, 
the system of nonlinear equations are directly solved using an 
iteration method, thus reduces the computational time involved in matrix 
inversions significantly. An assumed set of initial node voltages are 
iteratively updated according to an error criteria based on the 
assumed and the calculated node voltages. This method requires only 
one matrix inversion per time step. However, the convergence of the 
iterative process requires that the sum of the absolute values of the 
elements of any row or column vectors of the system Jacobian matrix be 

less than or equal to one .This conditon can be satisfied by choosing a 
proper time step. 


The method is applied to a harmonic converter circuit. This 
circuit consists of a sampler and a driver circuit. The sampler circuit 
uses nonlinear GaAs Schottky diodes as sampling switches. The driver 
circuit uses a step recovery diode for generating a narrow sharp 
triangular pulse to drive the sampler circuit. 

1.4.1 Time-domsin dei/ice modelling 


The characterization of a nonlinear device by a time-domain 
model is appropriate since the output parameters are related to input 
ones through device physics. Modelling of GaAs Schottky barrier diode 
and step recovery diode are considered in the present thesis. 


The modelling of GaAs Schottky diode has been reported by 
Thomas. W. Crowe ElOl. He used this model to study conversion loss in a 
mixer diode. Sussmon fort Eli] has developed a SPICE model for GaAs 
Schottky diode. A simple nonlinear model is used in this thesis. This 
model takes into account the series element and parasitic elements. 

The modelling of step recovery diode has been reported by 
S. Hamilton E123. He used a simple idealized model in harmonic generation 
circuits. The improved physical modelling of step recovery diode has 
been developed by John.L.Moll and Hamilton EiSl. The present model for 



^,,aru diode include the finite recombination effect,the finite 

step recove 

• ^ time and the series resistance. 

transition tif" 


1.5 


of the Thesis 


The major contribution of this thesis is an improved time- 
„,iusis method for nonlinear microwave circuit analysis, which 

domain a^aiy 

^p^^g,tionally more efficient and is specially designed for 

,-ircuits with high order nonlinearity. The method is also 
analysing ci . , ^ t. 

<)f handling signals with large number of harmonics . 


chapter 2 , 


the methods of analysis available in the open 


- i^nr the nonlinear microwave circuits are reviewed. A more 

literature 

^ /description of the proposed time-domain analysis method is 
detailed 

ented along with a discussion on the convergence criteria of 


the method- 


Schottky 


Hodellif^S' of nonlinear devices such as step recovery diode, 
diode and distributed elements such as lossy dispersive lines 


in chapter 3. It may be noted that a microwave circuit may 


is discussed 

/^rtlinear devices and linear elements interconnected bu 

contain 

line sections which may be dispersive. It is important to 

transmission 

incorporate all these types of elements in a microwave 


be able 


circuit analyse® 


software. 


An 


application example is presented in chapter 4 which 
use of nonlinear circuit analysis method in analysing a 


illustrates 

. ^nnverter circuit 
harmonic conv 


analysis method developed, enables one to evaluate the 
of a given circuit. However, if one wants to design a 

perf orrrianc? 

ji-h optimum performance, it is desirable to couple an 

circuit with 

orocedure with the analysis program, which automatically 

optimization P' 



adjust the circuit parameters for best performance. Such an attempt is 
made in chapter 5. Various aspects of this problem are discussed in 
some detail. 


In the concluding chapter, the work reported in this thesis is 
summarized. Some suggestions are made regarding the possible scope 
for further work. 
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2ANALYSIS OF GENERAL MICROWAVE NONLINEAR CIRCUITS 


2.1 Introduction 

Nonlinear circuit analysis methods can be classified as 
frequency-domain, time-domain or hybrid (mixed time and frequency- 
domain) methods depending on how the linear and nonlinear elements are 
analyzed. The linear elements including distributed components are 
easier to analyze in frequency-domain method , whereas the nonlinear 
elements are easier to simulate in time-domain. The hybrid methods have 
both these features but these methods require calculation of Fast 
Fourier Transforms and hence, these methods are limited to 
harmonically related signal components. A brief discussion of all these 
methods follows. 


2.2 Freouency-domain analysis 

Harmcxiic balance method : This method Cil is an extension of the 
Phaser analysis from linear to nonlinear differential equations. The 
solution of the nonlinear differential equations is assumed to be 
consisting of a linear combination of sinusoids, which are all 
harmonically related, making the solution periodic. When this assumed 
solution is substittfed into the system equations, the equations can be 
factored into a sum of purely sinusoidal terms, then superposition (due 
to the linearity of addition) and the orthogonality of sinusoids at 
different harmonics can be exploited to break up the resulting 
algebraic equations into a collection of simpler nonlinear equations, 
one for each harmonic. These equations are then solved by finding the 
coefficients of the sinusoids in the assumed solution that result in the 
balancing of the algebraic equations at each harmonic. 


A nonlinear differential equation has the form. 


f(xJcAJ)=0 


....( 2 . 1 ) 



where u e pCTo) is the stimulus waveform, x is the unknown waveform to 
be found and f is continous and real. p(To) denotes the set of all 
periodic funtions of bounded variation with period To, ie. u e p<To> 
implies that the function u satisfies the Dir ichlet- Jordan criterion [173. 
Assuming that the solution x exists, is real and belongs to p(To), then 

x(t) = 2) xdcte^*^*' , ....(2.2) 

K=-«’ 

where Uq = Zk/Tq. 


Substituting the assumed solution and its derivative into f, we 
get the resulting equation as Fourier series, 

f(x(t)Jc(OAJ(t)) = F(X,Ule)e’^^ ....(2.3) 

K=-“o 

where X = C....,X(-D,X(0),X(D,...J^ 

U = C....AX-D,U(0>,U(D,...J^ 

UCt) = ^ U(k)e*“°^. 

K=-» 

Finally the system of nonlinear equations have to be solved for 
X. 

F(X.U,k) = 0 for all k ....(2.4) 


To show how harmonic balance would be applied to circuit 
analysis, cnsider a circuit consisting of only voltage-controlled 
nonlinear resistors and capacitors, linear elements (lumped or 
distributed), and independent current sources. Assuming that the circuit 
has a periodic steady state solution and if the circuit has N nodes, 
then it can be described by 

t 

f(v,t) = i(v(t)> + J^cKv(l» + Jy(t-T)v(T)dT + is(t) = 0 ....(25) 

-OO 


where t is the time; v is the vector of node voltage waveforms; is is the 
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vector of source currents; ijq are differential functions representing, 
respectively, the sum of currents entering the nodes from the nonlinear 
resistors, and the sum of charges entering the nodes from the 
nonlinear capacitors; y is the itnpulse response matrix of the circuit 
with the nonlinear devices removed; and f is the function that maps the 
node voltage waveforms into the sum of the currents entering each 
node. 


It is necessary to transform both f(v) and v into the frequency- 
domain. To make the problem numerically tractable, the number of 
harmonics considered has to be some finite number. 

Since y is linear, the Laplace transform may be used to 
transform it into the frequency-domain y i — » Y. Also, since v is 
periodic and the circuit is stable 

t 

Jy(t-T>v(T>dT YV 

-OO 

where V contains the Fourier Coefficients of the voltage at each node 
and for each harmonic and Y is a block node admittance matrix for the 
linear portion of the circuit.. 


Y = CYmn^ m/i t Cl,2,..>0 

Ymn = CYmnOan W e 

where m/i are the node indices; k,l are the frequency indices. 

Ymnflc^ = f YnwiCikuo) if k = 1 

10 if k 1 

Now, equation (2.5) can be transformed into the frequency-domain. 


F(V) = I(V> + jftQ(V> + YV + Is = 0 


....( 2 . 6 ) 
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where Is is the Fourier Transform of is- Similarly, F<V), I(V), Q(V) are the 
Fourier Transforms of f(v), i(v), q(v> respectively. 

Cl = Ci^mrPi m/1 e {1,2,. ..,N} 

Anm = f ciiag{0/>)o,....,<H-l>a)o>> m = n 
I 0, m n. 


The nonlinear equations (2.6) have to be solved for each harmonic 
to determine the Fourier coefficients of the voltage at the each node. 

2.3 Hybrid (or mixed freauency and time-domain) method 

This method uses time-domain for nonlinear portion of the circuit 
and frequency-domain for linear portion of the circuit. A given circuit 
is partitioned into subcircuits, one subcircuit consisting of linear 
elements and the other one consisting of nonlinear elements only. The N 
branches of the linear-nonlinear interface connect the two subcircuits 
and define corresponding nodes. The current flowing out of the linear 
circuit must equal flowing into the other. Matching the frequency 
components in each branch satisfies the continuity equation for current 
(fig 2.1). 

.1 

The nonlinear circuit is represented by a set of nonlinear 
equations, 


ijTt) = gCVi(tj,....,v,gCt») ....C2.7) 

where g is an arbitrary nonlinear function and ij and vj are the Jth 
branch current and voltage respectively. 


The linear circuit is represented by an N by (N+M) matrix, 
obtained by standard linear circuit analysis. The M additional variables 
are the additional external nodes at which applied voltages (or 
currents) are present. The linear circuit matrix is calculated at each 
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Fig. 2.1 Analysis of Nonlinear Circuit by Hybrid Method. 




frequency component in the circuit. This matrix is 


V^CkW) 


HiiPtU) HiaCicU) 
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ilCWiJ) 

• 

— 

- 



* 


• 


Vn+i(“) 






. 

mm 


^ — 




1 

L- — 


....C2.8> 

where k = 0,1,2,... ,q and are the impedance or transfer variables, 

depending on which of the variables are voltages and which are 
currents. 


The initial guesses are established for the current phasors i^ku) 
at the interface branches at the dc, fundamental and harmonic 
frequencies (k = 0,1,2,.. The hybrid matrix for the linear circuit HOcu) 
is calculated at dc, the driving frequency u, and each harmonic. This is 
used with ijOcu) and the applied voltages in equation (2.8), to calculate 
the unknown phasor components of voltages at each of the N branches. 


Then using an expression, 


V j(t> = Re 


tc^o 


Vj(ku) e 


jkut 


....(2.9) 


to derive the time values of the branch voltages at times t = Ts,2Ts,..- 
and a similar expression for derivatives, the time samples of voltage 
and its derivatives are calculated at each of the N nodes. 


The values of ij(t) in the nonlinear side of the interface 
branches are obtained at corresponding time instants by substitution of 
the time samples of the voltage Vj(t) and its derivatives into eqn.(2.7) . 
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Using the DFT, the harmonic phasor components ij(ku) can be 
extracted. 

An error function is formed to compare the 'nonlinear' current 
estimates ij with the 'linear' estimates 


^l^lijCku) - i^<kti))p + ... |i^<k(t>> 


K=0 


...( 2 . 10 ) 


The error function is then minimized by forming new guesses for 
the current phasors ijOcu) from the old estimates, and repeating the 
above procedures until the error function lies below some threshold 
value (near zero). The fixed point method of Hicks and Khan C43 is used 
to force the error function to zero, by allowing the phasor currents to 
more closely approximate their true values on successive iterations. 
After rth iteration process, the current in the jth branch is 


i^(t) = ij^Ocu) e*“*^ 

ijp(t) = i^Ocu) 


The next iteration is then carried out with ij(p^jLj(kti>) formed by 


^jCr+D^“^ = pijpOcu) + (1-p) ijpOcu) 


...(2 .ID 


where p is determined by convergence considerations and 0<p^i. 


2.4 Time-domain av^ysis 

The time-domain method analyze the circuits at discrete time 
steps. The voltages and currents are described by a set of points over 
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ther period obtained at discrete time instants. The time-domain method 
involves the following steps. 

i) Converting the nonlinear differential equation into nonlinear 
algebraic equation by numerical integration of time equation. 

ii) Formulating the system matrix equations using standard 
circuit analysis such as nodal analysis, Tableau approach. 

iii) Solving the system equations by linearizing them .This 
nonlinear iteration process is repeated until convergence occurs. 

iv) The whole procedure is repeated for all time steps. 


Any time domain method uses the following algorithm: 

Algorithm: (Time incremental circuit analysis.) 

Input: Circuit description, Maxtime 
Output; Voltages, currents as functions of time. 

Start 

Nowtime = 0 
Repeat (T ime loop) 

Repeat (Nonlinear loop) 

Stamp circuit elements into matrix 'A' and 'S' 

Solve linear system AX = B. 

Until Nonlinear loop converged. 

Write (Nowtime, X(Nowtime)) 

Nowtime = Nowtime + At 
Until Nowtime = Maxtime 

End. 

The circuit equations can be formulated as a system of 
differential equations 

f(xpc,t) = 0 ...(2.12) 

where x is a vector of N circuit variables, and t is the time. 

|y ■ ' ■ ■ 
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The time-domain approach consists of discretizing the time 
interval into time points tj^, k = 1^,..., M and by replacing x at each time 
point by an integration formula. Usually backward differentiation 
formula is used. This approximate as a function of x at the 

previous time points. 


S ®P*k+l-p 

R 


...(2.13) 


where is the computed value of h = - tj, is the 

integration step size at time t|^ and P is the approximation order of the 
differentiation formula. The integration coefficient Op are chosen such 
that x^j = x(tjp, k= 1,2,. .J<, when the solution x(tj^ is a polynomial of 
degree P. For example, the Backward-Euler formula is obtained by 
choosing P=l, aQ=l, «£=1 for fixed time steps. By combining equations 
(2.12) and (2.13), the circuit equations can be expressed as 


which is a set of nonlinear algebraic equations. This new set of 
equations has to be solved at every time step to compute 


General purpose circuit analyzers like SPICE solve the nonlinear 
system of equations by the Newton-Raphson algorithm at each time step. 
This uses linearized equivalent circuit numerical models obtained from 
Newton-Raphson algorithm for nonlinear elements in each time step. The 
Newton-Raphson algorithm constructs a sequence of vectors *)£+l^ 
which eventually converges to Xfc+i- The vectors are the 

solutions of the set of linear equations; 

Bk.,1 <21® 

where the coefficient matrix and the vector ®k+l^ are related 

to the jth Newton-Raphson iteration at the time point At the each 

Newton-Raphson iteration process,the coefficient matrix' changes 



iG 

and hence the matrix A has to be inverted in each iteration ,to get 
^k+1^- several matrix inversion per time step can be avoided if 

we can solve the system of nonlinear equations directely.The direct 
nonlinear iteration process used in this thesis is explained in section 
2.4.2. 

2.4.1 Formulaiion of nodal admitt^vue matrix and eouu^alent current 
source i/ector 

The nodal analysis is used to formulate the system matrix from 
circuit equation as 

[ A ] [ V ] = [ B ] ..(2.16) 

where V is a set of node voltages, A is a node admittance matrix and E 
is a equivalent current source vector .Only linear resistors, linear 
capacitors,linear inductors, voltage and current sources and two 
terminal nonlinear elements are considered here. 

Linear resistor : If 'i' is the positive node and 'j' is the negative node 
of a resistor of value R,then 

A(ijl = A<i4l + 1/R 
Aaj) = Ati j) - 1/R 
A(j4) = - 1/R 

A(jj) = ms + 1/R ..<2.17) 

If one of the nodes, say 'i' is connected to a voltage source, 

then 


A<j4) = A(3 j) + 1/R 
B(j> = B(j) + Vs/R 


.42.BJ 


where Vg is the voltage source at the corresponding time instant. 


Linear capacitor : The Backward-Euler alforithm for solving the first 
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order differential equation v' = f(v) with a step size h is given by 

= vn + h f(v^j^) = Vp + h v'^m-I ..<2.19) 

where v* denotes differentiation with respect to time. For a linear 
capacitor, 


If we approximate the exact solutions v'Ctj^^^j,) and i<t^+'L) by and 

respectively, then 

..( 2 . 20 ) 

Substituting eqn.(2.20) in (2.19), we have 

= ^''n+i^ “ Cvn/h ..(221) 


This relationship can be represented by the equivalent linear 
port shown in fig. 2.2.(i) with port current i^+i port voltage 

Vpj^j^.The value of R remains fixed (assuming fixed step size), the current 
source depends on the voltage Vn; which is determined during the 
preceding time step. 

Linear inductor : The Backward-Euler numerical algorithm for solving 
first order differential equation i’=f(i) with a step size h is given by 

in+i = = in + ^’n+1 ...(222) 


For a linear inductor , v(t> = L i'(t) 
and this can be approximated by 

iVl = - <2^> 


Substituting EQn(2.23) in (2.22), we have, 



1 © 


i 



Fig.2.2(i) Discrete Circuit Model Associated 
with B-E Algorithm for a Linear 
Capacitor. 


in-fl 



Rg.2.2(ii) Discrete Circuit Model Associated 
with B-E Algorithm for a Linea r 
Inductor, 
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Wi = + in - <224) 

This eqn. can be represented by a linear equivalent one port 
shown in Fig. 2.2(ii). This circuit is the discrete circuit model 
associated with the B-E algorithm for a linear inductor. 

Current Sources - If the current is directed from node, i to j, then, 


BO) = BCD - Is 

BCj) = BCJ) + Ifi ...(225> 


where Is is the value of the current source. 

Nonlinear tuo terminal elements The nonlinear elements are 

characterized by 

i = f(v,v,t) ..(2J26) 

where v is the voltage across the device and the dot denotes 
differentiation with respect to time. This can be considered as a 
current source if the value of the voltage across the device is 
assumed in each iteration process. 

The formulated system of equations has the form 

[a][V]=[b] 

if a fixed time step is chosen, then, V can be written as 

[ V ] = [ A [ B ] ...(227) 

where A”^ is the constant matrix for a fixed time step and B a vector 
matrix function of nonlinear node voltages. This equation is then solved 
using direct iteration process, discussed in the following section. 
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2.42 Iteration method 

The set of nonlinear equations obtained from numerical 
integrations can be written as 

Vi = #i(Vi,V2,....Vn) 

V2= ^2<Vi,V2,....V„) 

Vn = ^n<Vi,V2'- 

or can be written as 

[ V ] = [ ^ ] ...( 2 J 2 ai 

where V is the node voltage vector and ^ is the nonlinear function of 
node voltages. 

A direct iteration method is used to find the actual root of the 
eQuations(2.28).The iteration process is 

^ ^ <n=0,i;2.....) ...C2J29> 

If this process of iteration converges ,then the limiting value is 
defininitely a root of the eq iation.(2.29).Let x=X be the solution of the 
equation and let = X + €„, where 6^ is the error in Xn.By expanding 
the right hand side of eqn.(2.29) in a Taylor's series, we have, 

fn+l “ •®ifn ' 

where ...(2.30) 

If a^^iO, then the errors resulting from the successive 
repetitions of the iteration process are re lated by 

fn+i = i €«+« " ®l*^n 
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The criteria for the convergence of the process is that laj = Vfixn 
must be less than one, and the magnitude of the error decreases 
exponentially with' n' increasing. 

We can make a better approximation to the initial guess in the 
process of iteration by 'exponential extrapolation' .If we consider only 
the first term in the Taylor's series expression(2 .3D),then 


or € 0^2 

If Xq is the initial guess and x^ and are the successive 
solutions, then 


(Xj - X) (Xq - X> = (x^- X)® 


therefore,an approximate root X* is given by 


X* = Xj - 


(Xg- 2Xt+ Xq) 


...(2.31) 


This process is repeated with starting Xp= X* till convergence is 
reached .An error function is formed to check for convergence .If x^ is 
the initial guess and x^ is the next iteration solution ,then error 
function is given by 


E = 




£ ***- X 
<■1 - 3 Cq 


0)2 



...(2.32) 


where. N is the number of nodes .The iteration process is continued till 
E lies below some threshold value .This limiting value for the iteration 
process is taken as O.iV, of the x value. 

Z.43Cor»^ergenceofi;heit6rationf>rocess 

The iteration process(2.29) converges if 
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or 


d#j(v) 

dVj 

j=t 


-A I d#.(v) 
jSL-f I dV; 

J=j ’ 


i: q;< 1 , tt = iA... .n) 


q;< 1 , (i = 1^2,. .../T) 


...(2.33) 


Consider the system of nonlinear equations 

A(h) V = Kv/O 

where A is a matrix, function of time step'h’,and B is a vector , function 
of node voltages and time step Therefore 'v' can be determined from the 
above equation as 

[V] =:[a]“^[B] = [C][B] ...(2B4) 

where C is the inverse of A. Now, v can be written in expanded form as. 









[V] = 



.+Cnnbn 

3 % 

^ntv,h) 


...(2.35) 


Consider a special case where,only and V2 are the node 
voltages For this,only bj^ and b2 are functions of and V2 

respectively,which are represented as 

bi = fi(Vi,V2) 

b2 = fzCV^.Vg) 


For a diode nonlinear device. 


= Is(exDaV-D + 


‘^jo dV 


(i-V/^)*" ^ 


b 2 = 


...( 2 .^) 
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where V = Vj^-V 2 - 

Applying the convergence criteria (2.33), 


8Vi 


3^4 


avr 


3b4 


ab4 


'ii ^ 


where 


ab 

a^ 


i _ 


= alsexpCaV) 




< 1 

uf C(z small' 
m(Vn-Vn_i>Cjo 


...(2.37) 


(i-V/#)"' ^ At (1-V/#)'^^ 



= alsexcKaV) 


")o 




m 


in(Vn-Vn-l>Cjo 
^ At (1-V/#)"*^^ 


The following approximations can be made to satisfy Eqn.(2.37) ; 

(i) Because of the exponential dc characteristic of the diode, the 
condition (2.37) is sometimes difficult to satisfy. To avoid this, a linear 
approximation to dc characteristic can be included at some voltage 
points and/or a piecewise linear model can be included. 

(ii) The junction capacitance can be taken as a fixed one, in which 
case, the capacitive element value goes to the A matrix. 

Finally, a proper time step is chosen experimentally to satisfy 


Eqn.(2.37). 
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3. MODELLING OF MICROWAVE DIODES AND TRANSMISSION LINES 

1 

3.1 Introduction 

The behaviour of a semiconductor device for a given geometry 
and doping profile is governed by the device physics^ and can be 
represented by a set of equations. The solution of these equations 
yields the device characteristic. The terminal characteristic of a 
device can be modelled by an equivalent circuit consisting of lumped 
components. The transmission lines have distributed characteristics. 
Modelling of transmission lines, only by lumped elements is not possible. 
But distributed components can be modelled in .time-domain by many 
sections of delay elements and lumped elements. 

32 Piecewise-Linear modelling 

Most electronic nonlinear device can be realistically modelled by 
a piecewise-linear characteristic. 

A nonlinear voltage controlled two terminal resistor is described 
by a canonical piecewise-linear represntation 1193. Consider an 
arbitrary continuous piecewise-linear function 'f' with 'p' distinct break- 
points (fig 3.1). Let m^,(i = 0,1,. ..,p) denote the slope of each segment. 
Vi<V 2 <.--<Vp. Then' f' can be described by a canonical form given by 

P 

i = f(v) = a + bv + Cjlv-Vj^l ...(3.1) 

i=l 

where the coefficients can be calculated explicitly from, 

b = (mjj + fHp) / 2 

C; = (m., - nv-j,), i = 1, 2,...,p 
P 

a = f(0> - ^ C; Nil- 
i=l 

3.3 Modelling of GaAs Schottky diode 


The nonlinear nature of a diode can be described by 




Fig .3.1 Piecewise-linear representation of a nwtlinear f motion 
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• i = f(v,v) ...(32) 

where v is the voltage across the diode and i is the current through 
the device and f is a nonlinear function describing the v-i 
characteristic . 


The dc current-voltage characteristic of a GaAs Schottky barrier 
diode is described by 

v/nvT 

i = Is C e -13 ...(3.3) 

where Is is saturation current, Vj = KT/q and n is the ideality factor, 
q- electronic charge, K - Boltzmann's constant and T - temperature. 


The saturation current is given by Richardson-Dushman equation 
E243 for the thermionic current. 

Is = AA*T^exp( ...(3.4) 

% 

where A - Junction area, A - Effective Richardson constant; T - 
temperature; - barrier height in volts. 

A* = (4ion*K^/h^ 

where m* = Effective electron mass, h - Planck's constant. 


The depletion or the junction capacitance of a Schottky barrier- 
diode is given by 


C. = 


Cjo 


(1-v/^)' 


.m 


...(3.5) 


where is the zero bias depletion capacitance, v is the applied 
reverse bias, > is the junction potential and m is the grading 



coefficient, m = i/2 for a step junction. 


3.3 .1 Eouh/alent circuit of SchoWof barrier diode 

The equivalent circuit for the GaAs Bchottky barrier diode is 
shown in fig. 3.2(i). Rs is the series resistance, Cp and Lg are the 
package reactances. The series resistance is a complex function of 
voltage and frequency and is determined usually by experiments. 


To avoid the occurence of overflow because of the exponential 
characteristic of the device, a linear characteristic tangent to the 
diode characteristic at some voltage point is included in the model. 
Similarly a linear characteristic is drawn at some point on the diode 
reverse bias characteristc to avoid underflow (fig 3.3). 


The dc current-voltage charcteristic of a typical GaAs Bchottky 
barrier diode is shown in fig 3.4. The parameters for this diode are 
l£=0.75 nA, n = 1.0, Vj = 26.3 mV, vJ = 0.2 V, CjQ= 1 pf and Rg = 12 fl. 
The c-v characteristic for this diode is shown in fig .3 .5. 

3.4 Modelling of step recovery diode 

The idealized equivalent circuit of the step recovery diode 
shows that the diode behaves as a two state capacitor of large 
capacitance under forward charge storage and small capacitance under 
reverse storage, with zero switching time between two states. The 
forward storage capacitance is the diffusion capacitance of the i layer 
and reverse capacitance is the depeltion capacitance of the space 
charge layer. 

A first order improvement in the diode idealized model is 
obtained by accounting for the recombination of the forward stored 
charge as expressed by the continuity equation. The recombination is a 
function of the effective minority carrier life-time Tj^. Bince forward 
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Fig .3 .4 i-v characteristic of a GaAs Schottky barrier diode 
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Fi3.3.5 c-v characteristic of a GaAs schottky barrier' diode 
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capacitance is determined by having a large number of free carriers 
there is a finite time constant t^ for switching from forward to 
reverse capacitance. 


In the frequency range (I/Tr) < f < (i/tj^), the idealized model is 
valid. At frequencies lower than i/T^, the finite recombination effects 
must be taken into account and at frequencies higher than i/t^^, the 
effects of finite switching time should be taken into account. 


The dc characteristic of a SRD is described by the relation 
shown in the eqn. (3.3). The saturation current is described by the 
relationship shown in eqn. (3.4). 


The forward diffusion capacitance is modelled as(Appendix-A) 




V/nVn 


nvn 


...(3.6) 


and the reverse depletion capacitance is described by 

C. 


Ct = 


JO 


(1-v/^)' 


.m 


...<3.7) 


3.4 .1 Equii/alent circuit of SRD 

The equivalent circuit for SRD is shown in fig .3 .2 (ii). Rs is the 
series resistance, Ls and Cp are the package elements. The overflow 
and the underflow problems are avoided by limiting the exponential 
characteristic at some point on the diode dc characteristic. 


The dc characteristic of a SRD is shown in fig .3 .6 and the c-v 
characteristic is shown in fig .3 .7, The parameters for this SRD are Is = 
O.i * 10“^® A, Vt= 26.3 mV, n = i, m = 0.5, vj = 0.6, = 0.5 pf, Tr = iO 

nS and Rs = 0.8 fl. 
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Fig.3.6 i-v characteristic of a step recovery diode 
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3.42 Charge storage and r^erse transient in SRD 

The large forward capacitance is due to the storage of electrons 
and holes in the central region of the diode. Not all of the injected 
charge can be removed because there is carrier recombination in the 
central layer during the storage period. Fig.3.6 shows the 'step 
recovery action' when the voltage is abruptly switched from forward to 
reverse bias. 


The total stored charge in the n-side of the diode is 


Qs 


qA J APndx 


...(3.8) 


where Wn is the n-layer thickness, A is the cross-section of the diode 
and Pn is the concentration distribution in the n-layer. 


The continuity equation for the hole current is 
1 ^ p - Po _ ap 

q sx Tp “ ■ at 


Multiplying the above equation by qA, we have 


A 


a^ 

dx 


+ OjA 


Ap 

Tp 


- qA 


§9 

at 




where Jp is the current density for holes, Tp is the life time for holes, 
Po is the concentration of holes at x=0 and p is the concentration at x. 


Integrating eqn.(3.9) from 0 to Wn and using eqn.(3.8) we have 

lp(0) - Ip(wn) = ^ ^ ...(3.10) 

where Ip is the diffusion current for holes and Qs is the stored 
charge. Eqn.O.iO) is called the charge control equation. Ip(wrP is a 
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Fig. 3. 8 Charge Storage and Reverse Transient in SR 

(i) The Driving Voltage Waveform. 

(ii) The Current and Charge Waveform. 

(ill) The Diode Voltage Waveform. 
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small current and can be neglected. 


The steady state forward current of the diode is obtained by 

setting = 0, 

ot 

If = Ip(0)=: ^ 


where Qsf is the stored charge under steady state forward bias 
condition. Thus 

~ If Tp ...C3.1.jD 


Now, when a reverse current is applied at t = 0 by reversing the 
bias voltage as shown in fig .(3 .8), the stored charge decays toward zero 
due to recombination and extraction. The charge control equation then 
becomes 


-Ip = 


at 


Tp 


...(3.12) 


where Ip is the reverse current. Solving this equation with eqn.C3.ii) as 
the initial condition, we get 

Qs<t) = TpC - Ip + + Ip) expC -t / Tp) 1 ...(3.13) 


The stored charge reaches zero at time tgCstorage time). 
Therefore using eqn.(3.13) we get te as 


ts = Tp In (1+I.p/Ip) 


...(3.14) 


The recovered charge is given by 


G|p — Ip tg 


.G.15) 


and the fraction of the total stored charge that is recovered is then. 




..(3.16) 



3B 


If the reverse current is much larger than the forward current, almost 
all of the stored charge is recovered. 


The injected charge to SRD is not completely uniform across the 
plane of the diode, so the 'fast transition' does not occur across the 
entire diode simultaneously. An empirical formula ti33 for this 
transition time is given by 

= E i - (79.5A/BRfo) 3^ is + VgR- 10~“/2.5 ...(3.17) 

where Vee is the breakdown voltage in volts, is is the delay time in sec 
and fo is the output frequency in GHz. This turn-on time is much faster 
than turn-off, so it has been neglected in the model. 

3.5 Modelling of lossy dispersit/e transmission lines 

The analysis and design of microwave circuits consisting of 
transmission lines can be facilitated by the use of equivalent circuits. 
The propagation characteristic, that is, the propagation constant and 
the impedence of ideal closed waveguides can be modelled in terms of a 
transmission line or a distributed parameter circuit consisting of LC 
elements. For ideal structures the series impedence and shunt 
admittance per unit length of the equivalent distributed parameter 
structures are rational functions and hence are easily realized in 
terms of LC elements. This, however, is not the case 
for open structures such as microstrips. Accurate evaluation of the 
propagation characteristic of these structures including dispersion, 
losses and impedance requires numerical techniques. The modelling of 
lossy dispersive transmission line is reported by V.K.Tripathi C203. He 
has developed equivalent circuits for lossy dispersive lines, which 
consist of ideal lumped linear elements and lossless TEM transmission 
lines. 


The equation describing non TEM lines are 
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dv 

3z 


az 


= L(f ) ^ + R i 

at 

= CW ^ + G V 


..X3.18) 


where L,C/? and G are elements per unit length of the line. 


The lossy dispersive transmission line and its equivalent circuit 
are shown in fig .3 .9 (a) and (b). In order to derive the expression for 
the series impedence, shunt admittance and the parameters of the ideal 
line, we can make use of ABCD parameters. 


(ABCD), 


'for the lossy JSf‘trSiv6 Ene ' 


cosh 71 21csinh71 
sii^ Cosh 71 

M3 


...(3.19) 


where Zq is the characteristic impedence, 7 is the complex propagation 
constant and 1 is the length of the line to be modelled. Zc and 7 are 
functions of frequency. 


(ABCDl 


for the e<i. hoaei ' 


i+ZY Z 

Y 1 


CosftzJnnn jZosinBoLm 


i Z*~ 



Y i+ZY 

I 

...(32D) 


where Zq, $0 and Lm are the impedance, phase constant and length of the 
ideal model line, and Z and Y are the series impedance and shunt 
admittances respectiv-'ely, of the two terminal lumped networks. 

The dispersive line is divided into a number of sections with 
length 1 whose maximum value depends on the frequency range of 
validity of the model and is chosen so that Z and Y are realizable in a 
desired form. Equating (3.19) and (3.20), [ since matrix elements A=D and 
AD-BC =11, we have 



(c) 

ig.3.9(a) Schematic of General Lossy Dispersive Line. 

(b) The Equivalent Circuit Model. 

(c) General Form of Lumped Circuits for the 
Series and Shunt Branches in (b). 




cosh 71 = (1+2ZY) CosiSoLm + i tZoY(l+ZY) 


= 2Y Cos^oLn, + j CI/Zq + ZoY^lsin^DUn 




The expressions for the freauency-dependent compls^ 
impedance and shunt admittance branch are derived from 
equations and are given by 


abo^® 


Zc 


- 1+ jsi 


cosh 71 - J 1+ jsinQSoLm) 


sinhTl 


sinhTl 


Y = 


sinhTl 


icos(j3oi-fn> - i i^i+isin(0oLm) Zq 


Zosin<ftol_m) 


..(3:2Z> 


the 


igngth 

The model parameters to be determined include tne 

. addition ^ 

phase constant, and the impedence of the ideal TEM line in 
the lumped circuits whose driving point impedance and admitt-en 


given by eqns.(3.22). 


The lengths of the line to be modelled, 1, and the 


n a desi*"®^ 

TEM line, Lm, are chosen such that Z and Y, are realizable i ks 

Q-f }l vC? 

form with ideal RLC elements. The maximum value of the line 
modelled, depends on the highest frequency at which the structu 
be modelled. 


p-ostef 


real CZU 


]f we choose to realize Z and Y in the general modifi®” 
form shown in fig.3.9(c>, then both functions must be positi'^® 
over the desired frequency band for which the model iS t-o ^ . tne 

f Kosized 

general form of impedance or admittance to be syntne=> 


modified Foster form is given by 
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Z or Y = 



(Cj + jci)) 


...( 323 ) 


where the constants B^, D| in the above equation determine the 

element values in the model. 
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4. APPUCATION OF THE TIME DOMAIN METHOD TO 
THE ANALYSIS OF HARMONIC CONVERTER CIRCUIT 


Harmonic converters are used in several microwave measurement 
systems such as network analysers, frequency counters, sampling 
oscilloscopes etc, because of their extremely broad band performance. 
The harmonic converter essentially consists of a sampling gate and a 
hold circuit, driven by a narrow sampling pulse. Obviously, in order to 
sample a microwave signal in the GHz frequency range, the drive pulse 
has to be necessarily much narrower than the period of the microwave 
signal. For the circuit to work as a down converter one has to choose 
the drive pulse repetition frequency such that 'n' times the frequency 
is equal to the signal frequency +IF or -IF frequency. By changing the 
harmonic number 'n', one can down convert any frequency in a multi 
octave microwave band with a small variation in the drive frequency. 
One major difference between the sample and hold circuit used in 
digitizing a signal at low frequencies and the harmonic converter is 
that the sampling rate is very much lower than that specified by 
Shannon's sampling theorem, but satisfies this criteria as far as the 
base band signal is concerned. 

4.1 Sampling Technique 

Sampling of a signal can be looked at as a multiplication process 
between the signal and the sampling pulses or the gate conductance 
function. 


y(0 = g(th x(t) ...(4.1) 

where x(t> is the microwave signal, g(t) is the conductance of the 
sampler gate and y(t) is the output of the sampling circuit. Fig. 4.1 (a) 
and (b) shows the sampling action. The gating signal has a wide 
frequency spectrum and hence, the sampling circuit has a wide 
frequency of operation. The sampling action in frequency domain can be 
modelled as a convolution. 



P‘g-4.1 (a) Sampling by Gating . 



Fig. 41 (b) Sampling Action 

x(t) is the microwave signal. 

g(0 is the conductance of the sampler gate. 

y(t) is the output of the sampling circuit. 
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Y(f) = G(f ) » X(f) 




where * denotes the convolution, given by 

oo 

G(f> * X<f) = J G(r) X(f-r) dr. 

— QQ 


In order to sample input signals in the GHz frequency range, the 
sampling pulse function has to be much narrower than the period of the 
highest frequency required to be sampled. This repetitive pulse train 
when transformed into frequency-domain gives a line spectrum which is 
extremely broad with each spectral line corresponding to a harmonic of 
the pulse repetition frequency. A harmonic converter can also be looked 
at as a mixer which generates an IF frequency taking the nth harmonic 
nearest to the signal frequency as LO frequency. But the analysis is 
not so simple. One needs to take into account the contribution due to 
the other frequency components as well. 

4.2 Sampler 

A harmonic converter circuit or a sampler circuit at microwave 
frequencies uses Schottky diodes as sampling switches, because of the 
fast switching characteristics required. The narrow sampling pulses 
are generated using a step recovery diode(SRD>. 


A schematic diagram of the sampler circuit is shown in fig .4 .2. D1 
and D2 are the two GaAs Schottky barrier sampling diodes. The sampling 
drive pulse is applied across the GaAs sampling diodes. Assuming no 
microwave input signal, the rising edge of the pulse voltage charges 
capacitors Ci and C2 to the peak pulse voltage through the sampling 
diodes D1 and D2. Then, as the drive pulse voltage returns to zero 
potential, the stored capacitor voltage appears across the diodes, 
reverse-biasing them. Resistor R3, R4 and R5 serve to discharge the 
capacitors slowly with the output voltage across Ri_ remaining zero all 
the time. Now, when a second drive pulse arrives, the GaAs diodes again 
become forward-biased, provided that the pulse peak voltage is greater 
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than the negative self-bias on the diodes. 


When a positive voltage is present at the microwave input to the 
sampler circuit, and the drive pulse is present, then the total current 
flowing through D2 is slightly more, and the total current flowing 
through Di is slightly less, than if no voltage were present at the 
microwave input. Therefore, capacitor C2 is charged slightly more and 
Ci is charged slightly less than if no voltage were present at the 
input. So, a signal voltage appears across the capacitors and causes 
currents to flow through resistors R3 and R4 and capacitors C3 and C4 
to load R^. If the input signal is negative, similiar operation causes the 
opposite-polarity signal to appear at the output. 


Ri and R2 simply terminate the microwave input transmission line 
in its characteristic impedance. GaAs diodes are chosen for minimum 
series resistance, junction capacitance, stray capacitance and 
inductance to have wide frequency of operation. 

4.3 Sampler driver . 

A SRD is used to generate a narrow sharp triangular pulse. The 
harmonic generator circuit used for this purpose is shown in fig .4 .3. 
This circuit when driven by a sinusoidal drive signal produces large 
negative-going triangular pulses of very short duration at the drive 
frequency. The SRD has a large capacitance under forward bias and a 
small capacitance under reverse bias. The diode stores the charge 
during forward bias and when the voltage is reversed, the stored charge 
starts discharging. When the charge goes to zero, it gives out a very 
sharp pulse at the output. 

In fig .4 .3, is the drive inductance and C-^ is the tuning 
capacitance. The section LI and Ci is used to match the network to the 
input transmission line characteristic impedance. Rg is the source 
resistance. In the output of the circuit, SRD is connected to a 


4 



I 


Fig. 43 Pulse Generating Circuit Based on a 
Step Recovery Diode . 
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microstrip transmission line. The narrow pulse generated by SRD circuit 
travels along the microstrip line which is coupled to a slot 
transmission line. The sampling pulse couples to the slotted line 
through the microstrip, that generates two opposite-polarity pulses to 
drive the GaAs Schottky barrier diodes (fig. 4.4). 


The SRD / Comb generator circuit produces a narrow pulse once 
per unit cycle, by storing energy in a drive inductance L^, just prior to 
the diode capacitance switching from forward to reverse capacitance. 
This energy appears across depletion capacitance C^, after switching, 
as a narrow pulse. The spectral content of the pulse can be found 
Fourier analysis. Approximating the pulse formed to a triangular pulse 
of height A, and width w and period T, we have 


f(t) = 0 ; 

= 2At/w + A ; 
= -2At/w + A ; 
= 0 ; 


-T/2 £ -w/2 

-w/2 i t i 0 

0 5 t £ w/2 
w/2 ^ t i T/2. 


The Fourier spectrum of this function is then, 

OO 

f (t> = 3^/2 + an cos(2itnl/T) , 

n=o 

where ao = Aw/T, = 2AT/n®*®w [l-cQS<rmw/T)3. The amplitude spectrum 
Cn is given by 

I Cn I = 2A/n*it®wf Cl-cos(rnrwf)l ...(43) 

The first zero occurs at the frequency = 2/w. 

The spectrum of sampler driver output for frequencies 900 MHz, 
950 MHz and w = 0.045 nsec are shown in fig. 4.5 and 4.6. These figures 
show that there is power in the harmonics upto 30 GHz. 
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Comb generator circuit design: The design of a comb generator circuit 
involves 1) Choosing appropriate SRD, 2) Drive inductance- capacitor 
Cy to tune out at fj^, 3) a matching network (LC section), 4) a bias 
resistor R^j. 


In the depletion or pulse formation interval, the output voltage 
across SRD is given by [121, 


eo<l) = 


Ii^niTc 

a|i - 


exp 




417 


sin jSt 


...(4.4) 


where Ii is inductor current at the instant of pulse formation. 

0 = J(i - / LoCo = 0)n 

€ = i/ 2 RL 41 ^ 

The pulse width tp is given by tp = it/i8 ...<4.5) 

Drive inductance can be found from tp as 

Lo = (tp/it)2/Co -..<4£) 

where Cp is the depletion capacitance of the diode. 

The tuning capacitor C,. resonates with Lp at therefore 

Lp = 1/bijCj 

From the above equation we have 

Ct = Co/(2fitp)* ...(4.7) 

The load resistor is obtained from 6 as 


Rl = tp/2ir€CD 


...(4.8) 


54 


For good output and pulse shape 6 is set at 0.3. 


The matching network is a two element network. The input 
impedance at terminals of Ct is a pure resistance which is given by 

Rin ~ 

If Rg is the source impedence and Rg /R^^ ^ then the values for the 
matching network are obtained from 


«-i = = ^<RgRin> /“i 

Cl = l/)Cci«i = i/[ J<RgRin) «i] ..X4.9) 

The important parameters of SRD are minority carrier life-time, 
zero bias depletion capacitance, cut-off frequency (f cut-off^ 
parasitic elements. 

= 1/2*RsCd -<4.10) 

where Rg is the series resistance .The diode with the minimum value of 
series resistance is chosen for the SRD circuit. 


4.4 An^ysis results 

4.4.1 SRD driller: 

The sampler driver circuit with designed values for input 
frequency 900 MHz and pulse width 0.045 nsec is shown in fig. 4.7. HP 
5082 - 0835 SRD is chosen for the circuit. The parameters for this SRD 
are Ig = 0.1 X 10'^® A, VT = 26.3 mV, n = 1, M = 0.5, Vj = 0.6, Cjo = 0.5 pf, 
Tp = 10 nsec and Rg = 0.8 Cl. To avoid overflow problems because of the 
exponential characteristic, the diode characteristic is linearized after 
0.8 V. Similarly underflow problem is avoided by limiting the actual 
characteristic to a linear one at -1 V. The series resistance is 
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included externally in the model for SRD. A time step of 0.1 nsec was 
chosen for the circuit. The calculated output waveform for many 
periods is shown in fig. 4.8. Similarly the SRD circuit has been solved 
for a input frequency of 950 MHz and pulse width 0.045 nsec. The 
calculated output waveform for this designed circuit is shown in fig. 
4.9. The input format for this circuit is given in Table 4.1. 

4.42 Sampler circuit 

The sampler circuit with element values is shown in fig. 4.10. A 
GaAs diode (HP 5082 - 0087) with parameters Is = 0.75 nA, VT = 26.3 mV, 
n = 1, M = 0.5, Vj = 0.2 V, Cjo = 1 pf and Rg = 12 ft has been chosen . To 
avoid overflow problems, the diode dc characteristic is linearized at 
0.3 V. The analysis of this circuit gave convergence problem. To avoid 
this, the junction capacitance of the diode has been taken as a fixed 
value of 1 pf. A time step of 0.005 nsec was chosen for this circuit. 
The circuit has been analysed upto 10 nsec. The input signal is varied 
from 1 to 24 GHz. Simultaneously, the SRD drive signal is also varied 
from 900 to 975 MHz so that the IF signal frequency is always fixed at 
100 MHz. The IF signal frquency is given by 
f IF = f sig - nf o 

where fgjg is the input signal frequency, fo is the drive frequency and 
n is the harmonic of drive signal. 

The computed IF output waveforms for different input 
frequencies 1, 2, 8, 24 GHz are shown in figs. 4.11, 4.12, 4.13 and 4.14 
respectively. The input format for this circuit is shown in Table 4.2. 
The conversion gain at different frequencies has been calculated and 
shown in fig. 4.15. 
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Fig .4 .8 Sampler driver cxjtput voltage waveform (900MHz ct'ive signal) 
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Fig.4.9 Sampler driver output voltage waveform (950MHz drive sigrral) 


100 

RCLVPDS 

10 6 10 1 0. lE-09 

V 1 0 

0.1E02 0.95E09 

R 1 2 

0.5E02 

R 2 0 

0.931E02 

C 2 0 

0.151E-10 

L 2 3 

0.185E-08 

C 3 0 

0.684E-10 

L 3 4 

0.41E-09 

R 4 0 

0-4775E02 i 

R 4 5 

0.8E00 
S 5 0 

O.lE-15 0.263E-01 O.lEOl 0.5E00 0.6E00 0.5E-12 0.1 


Table.4.1 Input ■Format for pulse generating circuit 
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Fig.4.ii Saiwler output voltage waveform for IGHz input signal 
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Fig.4.12 Sampler output voltage waveform for 2<Mz input signal 
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Fig .4 .14 Sampler fXJtput voltage waveform for 24GH2 input signal 
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RCLVPDS 

18 9 3 2 0 0.5E-11 

V 1 0 


0 . lE-01 

P 2 0 

O.IEIO 


0 . 1E02 

P 3 0 

0.9E09 

0.45E-10 

-0 . 1E02 

R 4 6 

0.1E03 

R 5 7 

0 . 1E03 

R 6 7 

0 . 22E05 

R 8 0 

0.1E05 

0.9E09 

0.45E-10 


C 2 4 

0 . 18E-11 
C 3 5 

0.18E-11 
C 6 8 

0 . lE-09 
C 7 8 

O.lE-09 
C 1 4 

0.3E-13 
C 1 5 

0.3E-X3 
C 1 4 

O.lE-11 
C 1 5 

O.lE-11 
D 4 1 

0.75E-09 0.263E-01 O.lEOl 0.5E00 0.2E00 0.2E-11 0 

D 1 5 

0.75E-09 0.263E-01 O.lEOl 0.5E00 0.2E00 0.2E-11 0 


Table 42 Input format for sanciler circuit 



Conversion efficiency.db 



Fig .4. 15 Conversion efficiency of a sampler circuit 


5. DESIGN OPTIMIZATICIN 


5.1 Introduction 

In any circuit or system design it is desirable that the circuit is 
not only functional but also optimum. The first task in optimization is 
to define what is optimum. This generally depends on the system and to 
some extent on the application. For example, we can define an optimum 
amplifier as the one that gives maximum gain while retaining the 
necessary bandwidth. It is well known that the Gain-Bandwidth product is 
constant in an amplifier, therefore maximum gain always occurs when the 
bandwidth is just sufficient or equal to that required by the constraint. 
Now, this can be a criteria of optimization in the case of an amplifier 
used as a power amplifier. But, if we are using the amplifier as the 
first stage of a receiver which contributes maximum to the degradation 
of signal-to-noise ratio, the optimality criteria is better chosen to be 
one of lowest noise figure subject to the bandwidth constraint. Thus, 
the criteria of optimization can depend on the application. 


Generally, two kinds of problems are encountered in an 
optimization process. One is the unconstrained and the other 
constrained optimization. The unconstrained problem is generally one of 
finding a maximum or minimum of a function. The procedure is to define 
an objective function which is minimized. (Maximization can be treated 
as minimizing the negative of the objective function.) 


The constrained optimization problem is normally solved by first 
converting the problem to an unconstrained one. This is done by 
constructing the objective function which includes the constraints 
automatically. There are several methods for doing this. The choice of 
the method depends on the problem and the type of constraints. 


Now, looking at the methods of minimization, one comes across a 
large number of algorithms, each one best suited for certain type of 
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problems or objective functions. The methods can be classified mainly 
into two categories, methods based on gradients and numerical search 
methods. Of course, the discussion here is confined to nonlinear 
function minimization. Linear functions are best minimized using the 
linear programming methods. 


All the gradient methods are based on the Taylor series 
expansion of the function about the minimum point. The drawbacks of 
these methods is the requirement of computation of gradients or the 
derivatives and their tendency to obtain local minimum or the minimum 
nearest to the starting point. There are some algorithms, such as 
Fletcher Powell algorithm, which do away with gradient computation but 
still suffer from the second drawback. 


When the objective function is well behaved smooth function or 
nearly a quadratic function, gradient methods perform very well, 
provided the initial point can be properly chosen. But when the 
objective function is complicated with large number of local minima, 
some of the numerical search techniques perform much better. Although, 
there is no sure sho . t way of determining global minimum the 
probability of hitting a global minimum can be sflldto be higher in some 
of the numerical search techniques .This, of course, is a qualitative 
statement supported only by experience. 


From these reasonings, a nonlinear minimization algorithm known 
as 'simplex method' developed by Nelder and Mead [331 was chosen for 
optimizing the harmonic converter circuit. A brief description of the 
method is given later in this chapter. 


As mentioned earlier in this section, the first task in finding an 
optimum solution or design, is to define or construct an objective 
function. All the literature on optimization deals with the problem of 
searching for the minimum, given an objective function ; and also deals 
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with how to convert a constrained problem to an unconstrained one. 
Therefore, the task of defining an optimum criteria and formulating an 
objective function is left to the individual, who is trying to design an 
optimum system, and is the best judge to decide on these matters, as he 
knows the system best. 


One of the major difficulties faced in the present problem of 
optimizing a harmonic converter circuit design is the computational time 
required. All the optimization algorithms are in general search for the 
minimum of the objective function iteratively, although search strategy 
is diffetxKvt in different methods. In order to reduce the computational 
time it is desirable that the objective function computation is fast and 
the number of function evaluations or iterations be minimum. The number 
of iterations required to find A minimum of the function depends on 
the minimization algorithm as well as the function behaviour with 
respect to the variables of optimization. The specific problem faced in 
optimizing harmonic converter circuit is one of inordinately large time 
for function computation. 


An important parameter of a harmonic converter is its 
conversion loss. It was sought to minimize the conversion loss by 
appropriate choice of component values. Therefore, the objective 
function would logically be the conversion loss and the circuit elements 
as the variabes of optimization. Now, to express conversion loss in 
terms of element values is not simple; in fact, complete time-domain 
analysis as developed in earlier chapters has to be gone through for 
several time steps to construct the output signal, in terms of which 
the conversion loss is defined. This takes large computational time. 


In order to reduce the computational time, it is desirable to 
reformulate the objective function. What is attempted in this thesis is 
to define a new function which in some way related to the conversion 
loss; and it is hoped that by optimizing this new parameter, which takes 
much less time to compute, we are actually minimizing the conversion 
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loss. There is no proof provided to show that this really happens. 
Instead of computing the output waveform once several time steps and 
then finding the maximum or the rms value of that waveform, a single 
sample value at a time step just after the sample pulse amplitude goes 
to zero is used in the objective function. The ratio of the input voltage 
at the sampling instant to the sampled output voltage value is taken 
to be some representative of the conversion loss. By heuristic 
reasoning one can see that this ratio in some way directly related to 
the conversion loss and minimizing this ratio, we are actually minimizing 
the conversion loss in some indirect way. 


52. ObjectjMe function and i/^iables of optimization for sampler circuit. 

The function of the sampler is to convert microwave frequencies 
to IF frequencies ;,withv. minimum signal loss. So, objective function for 
sampler circuit is the conversion loss. The objective function is 
defined as 


F ...(5.1) 

where Zq is the characteristic impedance of the input line and l^is the 
load resistance. 


The ideal diode for sampler circuit is a diode with zero series 
resistance, parasitic elements and zero cut-off potential. So, a GaAs 
diode with minimum series resistance, parasitic elements has been 
chosen. The parameters for the variables of optimization process taken 
are R3, R4 and R5. The SRD drive input with minimum pulse width gives 
flatter frequency response, but greater conversion loss. 


The determination of conversion loss function in the analysis 
involves a nonlinear problem. The behaviour of the objective function is 
totally unknown and is highly nonlinear. It can have multiple minima. 
Such functions don't easily yield to gradient methods because, these 
methods generally give the minimum nearest to the starting point, near 


the global minimum which is often difficult in the case of nonlinear 
functions . 

5.3 Numerical search algor ithm for function minimization 

The geometrical figure formed by a set of (n+i) points in an n- 
dimensional space is called a simplex. The 'n' points forms a set of 
variables of optimization. 


The objective function formulation for conversion loss is 
nonlinear and can have multiple minima. The simplex search method of 
Nelder and Mead [333 which is a powerful method in handling nonlinear 
functions has been used in this thesis. This method has multiple 
starting points and has greater probability of locating the global 
minimum. This starts with a set of (n+i) points, termed simplex, the next 
iterative point being formed by reflecting the point with the highest 
objective function value about the centroid of the rest of the n-points. 
The line joining the centroid and the point with the highest function 
value gives the general direction in which the function decreases and 
the simplex is expanded or contracted in the direction depending on the 
function value of the new calculated point . The search is terminated 
when the difference between the highest and the lowest function values, 
among the points of the simplex, is less than a pre-specified value. 

5.4 Computational aspects of minimization 

The evaluation of the objective function formulated as in eqn.(5.i) 
takes large time, that is, around 38 minutes in PC/XT, for a single 
function evaluation. So it is not practically possible to implement this. 
Therefore a time point after the SRD drive signal falls to zero in the 
first cycle is assumed to represnt the maximum value of the output, IF 
signal. The objective function is then defined as 

F 0 05 ns)/Zo 

" 0-05 ns)/R,_ 





Calculate initial values Pj^ and 
■Determine Calculate p h) 

Determine pj^ and yj^ 

» _ 

Form p = (l+a)p - ap|^ 

Calculate y* 



Replace P|^ by p 
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which is a representative of the conversion loss in some indirect way. 

The variables for the optimization taken are R3, R4 and R5. The 
search simplex is formed around the point (R3 = iODft, R4 = iOOA, R5 = 
22 K) by symmetrically perturbing each component by a 'fixed amount to 
give a new point. In the process of function mininizatioa the variables 
R3, R4 and R5 may go to negative value. Whenever the variable is going 
to negative value, the variable is changed to some positive value and 
the process is continued (R3, R4 and R5 are taken as iOfi, iOA, 22K 
respectively). The minimization is terminated when ty^ ■ Wil < 0.001 is 
obtained. The process took 26 iterations for minimization. The optimized 
variables are R3 = 12.4 A, R4 = 12.4 A and R5 = 27.95 K for 1 GHz input 
signal. The objective function value in dB for the optimized sampler 
circuit is 4.55 dE for 1 GHz input signal. 


The elements Cl and C2 are also taken as variables and tried in 
the simplex program, but the objective function behaved randomly and we 
could not get the minima for this case. 

5.5 Conclusion and discussions 

The iteration time in the optimization process is found to depend 
on the constraint level .The number of iterations also depend on the 
number of variables. 


The Schottky diode with minimum series resistance and parasitic 
elements was chosen .The ideal diode for the sampler circuit can be 
easily determined and so the diode parameters were not included as 
optimization variables .Also,f or SRD driver circuit,the output voltage has 
to be a delta function for generating a narrow sharp triangulur pulse 
and it can be obtained if the depletion capacitance of the diode is too 
small and the forward diffusion capacitance is high .Therefore the 
pulse generating circuit was not optimized . 
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6. CXJNCLUSION AND SCOPE FOR FURTHER WORK 

6.1 Present work 

In this thesis, a time-domain method for nonlinear microwave 
circuits has been suggested and successfully applied to harmonic 
converter circuit. The method is efficient for solving a class of 
strongly nonlinear microwave circuits and/or input signals having a 
large number of harmonics. The system matrix is formulated using nodal 
analysis and the matrix diagonal elements are all strong. The input 
source variables are avoided by the present formulation. The method 
requires only one matrix inversion per time step. 


The modelling of Schottky diode and step recovery diode have 
been discussed. Also the modelling of lossy dispersive transmission line 
is discussed. The harmonic converter circuit has been solved by the 
present method and results have been obtained for different 
frequencies. Finally, a possible application of function minimization 
algorithm to optimize the performance of the circuit has been 
discussed. 

SZ Limitations of the present method and scope for further work 

The present time-domain method requires a convergence criteria 
as discussed in the section 2.4.3. This condition, sometimes, is difficult 
to satisfy. So, present method fails for those nonlinear circuits, for 
which convergence criteria cannot be satisfied. Under these 
circumstances, one or more Newton-Raphson iterations can be applied, 
whenever the present iteration process diverges. 


The transmission line (lossy and dispersive) model has not been 
included in the present work. This model can be included and a wide 
class of microwave circuits can be analyzed. Also, models for other 
nonlinear devices such as PIN diode. Varactor, Tunnel diode, Gunn diode 
and MESFET can be included in the program. 



APPENDDC-A 
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Device Physics 

1. GaAs Schottky barrier diode 

i) Depletion or junction capacitance: The space charge layer 
width of a Schottky barrier on a uniformly doped n-type GaAs 
semiconductor is 




2 t<^-v) 




where £ is the permittivity of the GaAs semiconductor, is the doping 
density, ^ is the contact potential and v is the reverse bias voltage. 
The capacitance of the junction is therefore given by 


c - ^ 

3 - Xj, “ M2(<f-v) 


A 


where A is the area of the diode. The above equation can be written as 


c- = -Ei^ 

^ Al(i-V/^) 


where = 


OeNr 


\ 2* 


A. 


2. Step recovery diode 

D Diffusion or forward capacitance : In the forward biased 
condition, the charge in the p-region of SRD is given by 


Qp = oALpPno Cexp(v/nV-r^ - D 
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and the charge in the n-region is given by 
Qn = - tlALnOpo (expCv/nV^.) - 1) 

where q is the electronic charge, A is the area of the diode , Lp and Ln 
are the diffusion lengths for holes and electrons, Pno and npo are the 
minority hole and electron concentration in n and p-regions 
respectively, V^- = KT/q and n is the ideality factor. The total diffusion 
charge is then given by 


Qhi-p-p ~ Qp I On I = ~ dAQ-nnpQ+LpPpjo) (exp(v/nVj) - 1) 

Therefore diffusion capacitance is 

^ ^ tV«Lnnpo+lpPno) exp(v/nVT) 

^ = ‘«diff/^ = riJZ 


Pno » npo. 


Cjj w 


qALpPno exp{v/nV-,.> 

hVt 


or this can be written as 

Cp = dsTp/nV-f) exp(v/nVT) 

2 

where Is = qADpPrK/I-p and Tp = Lp /Dp. Tp is the life time for holes in 
the n-region and Dp is the diffusion constant for holes. 


ii) Depletion or transition capacitance : The charge in the space 
charge layer during reverse bias condition is given by 


Qj = A 




Zc-gjCj' 


N^+Nu 


- v) 
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where A is area of the junctioafig; is the permittivity of silicon,^ is the 
contact potential, is the acceptor concentration in p-region,Nj) is the 
donor concentration in n-region and v is the applied reverse bias 


The transition/depletion capacitance is therefore given by 


d^i 






Z <^ - v) 


or can be written as 


^30 

^ a1(1-v/<^) 


where Cj is the zero bias capacitance. 
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